EXISTENCE AND STABILITY OF STANDING PULSES IN NEURAL 
NETWORKS: II. STABILITY 



YIXIN GUO* AND CARSON C. CHOW+ 

Abstract. We analyze the stability of standing pulse solutions of a neural network integro- 
diffcrcntial equation. The network consists of a coarse-grained layer of neurons synaptically connected 
by lateral inhibition with a non-saturating nonlinear gain function. When two standing single-pulse 
solutions coexist, the small pulse is unstable, and the large pulse is stable. The large single-pulse 
is bistable with the "ail-off" state. This bistable localized activity may have strong implications 
for the mechanism underlying working memory. We show that dimple pulses have similar stability 
properties to large pulses but double pulses are unstable. 
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1. Introduction. In the accompanying paper [25], we considered stationary lo- 
calized self-sustaining solutions of an integro-diffcrcntial neural network or neural field 
equation. The pulses are bistable with an inactive neural state and could be the un- 
derlying mechanism of persistent neuronal activity responsible for working memory. 
However, in order to serve as a memory, these states must be stable to perturbations. 
Here we compute the linear stability of stationary pulse states. 

The neural field equation has the form 

(1.1) + u{Xj t ) = |°° W (x - y)f[u(y)]dy 
with a nonsaturating gain function 

(1.2) f[u] = [a{u{y, t) - u T ) + l]Q(u - u T ) 

where 0(-) is the Heavisidc function, and "wizard hat" connection function 

(1.3) w(x) = Ae- a W-e^ xl . 

In Ref |25| . we considered stationary solutions uq(x), where uq(x) > u T on an 
interval —x T < x < x T , u(x T ,t) = u T , and u(x,t) = uq(x) satisfies the stationary 
integral equation 

w(x - y)[a(u (y) - u T ) + l]dy, 

-XT 

We have shown the existence of pulse solutions of equation l|1.4|) in the form of single- 
pulses, dimple pulses and double pulses [^11221 ■ Examples can be seen in Figs. II. ll 
17.11 and 16.11 We constructed the pulses by converting the integral equation 1)1. 4f) into 
piecewise-linear ODEs and then matching their solutions at the threshold points x T - 
When the excitation A and gain a is small, there are no pulse solutions. If cither 
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Fig. 1.1. Single-pulse solution. 



is increased, there is a saddle-node bifurcation where two coexisting single-pulses, a 
small one and a large one, arise. As the gain or excitation increases, more than two 
pulses can coexist. For certain parameters, the large pulse can become a dimple pulse, 
and a dimple pulse can become a double pulse OJ • 

In this paper, we derive an eigenvalue equation to analyze the stability of the pulse 
solutions. While our eigenvalue equation is valid for any continuous and integrable 
connection function w(x), we explicitly compute the eigenvalues for l|1.3|l . For the 
cases that we tested, we find that the small pulse is unstable and the large pulse is 
stable. If there is a third (larger) pulse then it is unstable. The stability properties of 
dimple pulses are the same as corresponding large pulses. Double pulses are unstable. 

2. Eigenvalue equation for stability. We consider small perturbations around 
a stationary pulse solution by substituting u(x, t) = uq(x) + ev(x, t) into where 
e > is small. Since the pulse solutions are localized in space, we must assume the 
perturbation to the pulse will lead to time dependent changes to the boundaries of the 
stationary pulse (i.e where uq(x t ) — u T ). Thus the boundaries —x T and x T become 
time dependent functions 



where eAi(t) and eA2(t) are the changes of the boundaries — x T and x T produced by 
the small perturbations. Inserting u(x, t) = uq(x) + ev(x, t) into (|l.lfl and eliminating 
the stationary solution with l|1.4|) gives 



(2.1) 
(2.2) 



xi(t) 
x 2 (t) 



-x T + eAx(t) 
x T + eA 2 (t) 



(2.3) 




where 



(2.4) 




-(i T +eAi(t)) 
/■XT+eA 2 (t) 




u T ) + l)dy, 



(2.5) 





u T ) + l]dy. 
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Expanding the integrals I\ and I2 to order e yields the linearized dynamics for the 
perturbations v(x, t) 

(2.6) vt(x, t) + v(x, t) = a / w(x — y)v{y 1 t)dy — w{x + x T ) Ai + w(x — x T )A2- 



The time dependence of Ai and A2 is found by using the fact that u(x,t) is equal 
to the threshold u T at the boundaries of the pulse. Inserting ((2.1(1 and 1(2.2(1 into the 
boundary condition u(x\(t),t) = u T and expanding to first order in e leads to 

(2.7) Ai(t) = - <~ Xt ^ 



where 
(2.8) 

Similarly, 



duo(x) 



dx 



> 0. 

X — — Xt 



(2.9) A 2 (*) = ^^. 

c 

Consider time variations of v(x,t) that obey 

(2.10) v(x,t) = v(x)e xt 

where v(x) is a bounded and continuous function that decays to exponentially as 
x -> ±00. Substitute (|2~TU|) with lf2~7ll and l(2~3jl into lf2~B|) . to obtain 

(2.11) (1 + X)v(x) = w(x - x t ) V ^ Xt ^ + w(x + x T ) V(y Xt ^ + a [ w(x - y)v{y)dy. 

c c J_ XT 

where A is an eigenvalue with corresponding eigenfunction v{x). Equation (|2.11|l is an 
eigenvalue problem that governs the stability of small perturbations to pulse solutions 
of the neural field equation ^l.ip . If the real parts of all the eigenvalues are negative, 
the stationary pulse solution uq{x) is stable. If the real part of one of the eigenvalues 
is positive, uq(x) is unstable. 

We define an operator L: C [—x T , x T ] — > C [—x T , x T \: 

v(x T ) v(—x T ) f XT 

(2.12) Lv(x) = w(x — x T ) \- w(x + x T ) ha / w(x — y)v(y)dy. 

c c J_ XT 

Then the eigenvalue equation H2.11JI becomes 

(2.13) (1 + X)v(x) = L{v{x)) on C[-x T ,x T ]. 

We show in the Appendix ( Theorem 18.7(1 that L is a compact operator. We also 
show the following properties of the eigenvalue equation p.lljl : 

1. Eigenvalues A are always real (Theorem 18.4(1 . 

2fco 

2. Eigenvalues A are bounded by Ab = h 2akix T ~ 1 where fco is the max- 

c 

imum of |w(x)| on [0, 2x T ] and \w(x — y)\ < ki for all (x.y) G J x J, 
J = [— x T ,x T ] fTheorem 18.5(1 . 
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3. Zero is always an eigenvalue ( Theorem 18 .6f) . 

4. A = — 1 is the only possible accumulation point of the eigenvalues (Theorem 
18.81) . Thus, the only possible essential spectrum of operator L is located at 
A = — 1 implying that the discrete spectrum of L (i.e. eigenvalues of l|2.11|0 
captures all of the stability properties. 

We use these properties to compute the discrete eigenvalues to determine stability of 
the pulse solutions. 

3. Linear stability analysis of the Amari case (a = 0). Amari 3 computed 
the stability of pulse solutions to for a = 0. He obtained stability by comput- 
ing the dynamics of the pulse boundary points. He found that the small pulse is 
always unstable and the large pulse is always stable. Pinto and Ermentrout 03] later 
confirmed Amari's results by deriving an eigenvalue problem for small perturbations. 

We consider a stationary pulse solution of Hl.lfl with width x T . Applying eigen- 
value equation I2.11fl to the Amari case yields 

(3.1) (1 + X)v(x) = w(x - x T )^l + w(x + x t ) v( -~ Xt ^ ee Ti(v(x)), 

c c 

where T\ is a compact operator on C[—x T ,x T ] (see Theorem 18. 7|) . The spectrum of a 
compact operator is a countable set with no accumulation point different from zero. 
Therefore, the only possible location of the essential spectrum for T\ is at A = —1. 
This implies that instability of a pulse is indicated by the existence of a positive 
discrete eigenvalue. 

The eigenvalue A can be obtained by setting x — —x T and x — x T in i|3.1|) to give 
a two dimensional system 

(3.2) (l + A - ^) v(x T ) ^-v(-x T ) = 

(3.3) -^-v(x T ) + (l + X -^B) v{ - Xt) = o 



This is identical to the eigenvalue equation of Ref 031- Setting the determinant of 
system H3.2(l and l|3.3(l to zero gives the eigenvalues 

(3.4) x ^ W (0)± W (2x T )_ h 

c 

which agrees with Rcf. |43|. 

The stationary solution of the Amari problem satisfies 

/XT PX — XT 

w(x - y)dy = j w(y)dy 
-XT J X-\-Xt 

Differentiating u{x) yields u'(x) = w(x + x T ) — w(x — x T ) implying 

(3.6) u'{-x T ) = w(0) - w(2x T ) = c. 
Inserting into l|3.4|) gives the eigenvalues 

(3.7) A = W(0)+ " (2 ^1, 
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The zero eigenvalue was expected from translational symmetry. Since w(0) > w(2xt), 
the sign of c alone determines stability of the pulse. Recall that the small and large 
pulse arise from a saddle node bifurcation At the saddle node bifurcation, 

both eigenvalues are zero. Thus, setting A = in (|3.7|) shows that the width of the 
pulse satisfies w(2x T ) = For our connection function, w(x) has only one zero 
at xo for w(x) on (0,oo) (see [21 EH). Thus x T = xo/2 at the saddle node. For the 
large pulse, x T > xo/2, implying w{2x T ) < and c > 0. Conversely, c < for the 
small pulse. Thus the large pulse is stable and small pulse is unstable. 

Consider the example: a = 2.4, A = 2.8, u T = 0.400273, a = 0. There exist 
two single-pulses, the large pulse 1 and the small pulse s |24l I25| . For the pulse 1, 
x l T = 0.607255 gives the non-zero eigenvalue A = —0.165986 < 0, indicating it is 
stable. For the small pulse s, x T = 0.21325, gives A = 0.488339 > indicating it is 
unstable. 

4. Computing the eigenvalues. For the case of a > 0, we must compute the 
eigenvalues of (|2.11(l with the integral operator. Our strategy is to reduce the integral 
equation to a piecewise linear ODE on three separate regions. The discrete spectrum 
can then be obtained from the zeros of the determinant of a linear system based on 
the matching conditions between the regions. This approach is similar to the Evans 
function method E2 El E3 d] . 

4.1. ODE form of the eigenvalue problem. We transform H2.11JI (with the 
connection function defined by into three piecewise linear ordinary differential 

equations (ODEs) on (— oo,x Tl ), (— x Tl ,x Tl ) and (— x Tll oo). The ODEs then obey a 
set of matching conditions at x = x Tl and x — —x Tl . 

On the domain x € (— x Tl , x Tl ), we can write (|2.11(l in the form 

(4.1) (1 + X)v(x) = T^x) +h-h + h-h 

where 



X 



h(x)=a Ae- a ^- y) v(y)dy, h{x) = a j e~ ( ~ x -^v(y)dy 

J -X Tl J— XT t 



X 



and 

(4.2) TAx) = w(x - x T , + w { x + x ) v (- Xt i\ 

c c 

Differentiating l|4.1|l repeatedly gives 

(4.3) (1 + X)v'(x) = T[[x) - ah +h + ah - h 

(4.4) (1 + X)v"{x) = T{'(x) + a 2 h - I 2 + a 2 I 3 - h + 2a{l - aA)v{x) 

(4.5) (1 + X)v"'{x) = T{"(x) - a 3 h + I 2 + a 3 I 3 - 7 4 + 2a(l - aA)v'(x) 

(4.6) (1 + X)v""(x) = Ti"(x) + a 4 h -I 2 + a 4 / 3 - 7 4 + 2a(l - a 3 A)v(x) + 

2a(l ~ aA)v"{x) 



where we have used 



I[ = —all + aAv(x), I' 2 = —h + av(x) 

I 3 = a/3 — aAv(x), I' A — I4 — av{x). 
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Taking l|4.5|l — a 2 (14. 1|) and rearranging gives 

(4.7) I 2 + h = -5- — [(A + l)v" + (2aaA -2a- a 2 A - a 2 )v + a 2 T x - T"l 
a 2 — 1 

Substituting (|4.7|l back into (|4.1() leads to 
1 



(4.8) 



h 



[(A + l)v" + (2aaA - 2a - X - l)v + Ti - T['] 



Substituting both 1)4. 7fl and 1)4. 8|) into l|4.6|) . results in a fourth order ordinary differ- 
ential equation for v on the domain x € (— ir r , a; T ) 



(4.9) 



1 + A 



(1 + A)(a 2 + 1) 



III! 



+ 2(1 - <zA) 



Using Ti"(x) - (1 + a 2 )Tf (a;) + a 2 Ti(a;) = (obtained by differentiating T x (x)) and 
simplifying, leads to 



(4.10) 



(1 + X)v -Bv" + Cv = 0, xe(-x T ,x T ) 



where B = (1 + A)(a 2 + 1) + 2a(l - aA), and C = (A + l)a 2 - 2aa(A - a). 
On the domain x G (a; r ,oo), (|2.11() can be written as 



(4.11) 
where 



(l + X)v = T 1 + J x - J 2 , 



Ji = aA 



-a(x-y) 



v{y)dy, 



J-2 



-(x-y) 



v{y)dy, 



and Ti is defined by (|4.2|l on the domain [x T ,oo). 

Differentiating H4.11|) and using J[ — —aJi and J' 2 — — Ji gives 



(4.12) 
(4.13) 



(1 + X)v' (x) = T{ - aJi + J 2 , 
{l + X)v"{x)=T{' + a 2 J 1 -J 2 



Taking a^T^ + (a + D KWi + ffTJ^i and using T" + (1 + a)T{ + aTi = leads to 

(4.14) v" + (a + l)v' + av = 0, x £ (x T ,oo). 
Similarly, the ODE on (— oo, —x T ) is given by 

(4.15) v" - (a + l)v' + av = 0, x e (-oo, -x T ). 

In summary, the eigenvalue problem i|2.11[) reduces to three ODEs: 

(ODE I) v" - (a+l)v' + av = 0, xe(-oo,-x T ), 

(ODE II) (1 + X)v"" - Bv" + Cv = 0, xe(-x T ,x T ), 

(ODE III) v" + (a + l)v' + av = 0, x e {x Tl oo) 7 



where B = (1 + A)(a 2 + 1) + 2a(l - a A) and C = (A + l)a 2 - 2aa(A - a). 
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4.2. Matching Conditions. The solutions of ODE I, II, and III and their first 
three derivatives must satisfy a set of matching conditions across the boundary points 
— x T and x T . We derive these conditions from the original eigenvalue equation l|2.11|l 
which we write as 

(4.16) c(l + X)v(x) = w(x — x T )v(x T ) + w(x + x T )v{—x T ) + caW(x), 



From Ij4.lt) [I . we see that v(x) 



w(x — y)v(y)dy, x € (—00,00). 

-XT 

is continuous on (—00,00). However, w(x) has a cusp at x — which will lead to 
discontinuities in the derivatives of v(x) across the boundary points — x T and x T . 

We first probe the discontinuities of W(x) and its derivatives. W{x) is continuous 
on (—00,00). By a change of variables, W(x) — f x ^ x T w ( z ) v ( x ~ z )dz, from which 
we obtain 

pX + Xt 

W'(x) — w(x + x T )v(— x T ) — w(x — x T )v(x T ) + / w(z)v'(x — z)dz, 



x — xt 



indicating that W'(x) is also continuous on (—00,00). However W'(x) is not smooth 
at — x T and x T . Differentiating W'(x) for x 7^ —x T ,x T gives 

W"(x) = w'(x + x T )v(—x T ) — w'(x — x T )v(x T ) + w(x + Xt)v'(-x^) 

r-X — XT 

— w(x ~ x T )v'(x^ ) — / w{z)v" (x — z)dz 

J x+xt 

where v'(-x^) — lim v'(x) for x > —x T (right limit) and v'(x~) = lim v'(x) for 

x < x T (left limit). 

Using the following convention: 



to represent the jump at the boundaries, we find that 

[W"(x T )\ = W"{x)\ x=x+t - W"(x)\ x=x - = - [w'(0)]v(x T ) 
[W"{-x t )] = W"{x)\ x= _ x+t - W"{x)\ x _ x - = [w'(0)}v(-x T ) 

We differentiate W"(x) for x 7^ — x T , x T and find 

[W'"{x T )] = - K(0)] v(x T )- K(0)] v'(x-) 
[W'"(-x T )] = K(0)] v(-x T ) + K(0)] v'(-xt) 

To find the matching conditions for the derivatives of v(x), we differentiate l|4.16|l 
with respect to x for x ^ —x T ,x T , and obtain 

c(l + X)v (x) = w (x — x T )v(x T ) + W (x + x T )v{— x T ) + caW'(x). 

v'(x) is discontinuous at the boundaries because of the discontinuity of w'(x) at x = 0. 
Therefore 

W(?t)]= h \,A w'(0)}v(x T ), 
c(l + A) 

Iv'(-Xt)] = ^-tt[w'(0)}v(-Xt). 
c(l + A) 
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Differentiating 14. f 6|) twice yields 

c(f + X)v"(x) — w"(x — x T )v{x T ) + w"(x + x T )v{—x T ) + caW"{x) x ^ —x T ,x T 

There are discontinuities of v"{x) at —x T and x T that come from W"(— x T ) and 
W"(x T ). Note that w"(0~) = w"(0 + ). The jump conditions of v"(x) at — x T and x T 
are 

W'[x T )] = ^ [W"(x T )] = -j^ K(0)] v(x T ), 
[v"(-x T )\ = ^ [W"{-x T )\ = ^ K(0)] v(-x T ). 
By differentiating a third time we find the jump conditions for v"'(x) at —x T and x T : 
KM] - ~7T~~< \T K"(0)] + ^ [W"'(x T )] 



c(l + 


A) 


f 




c(l + 


A) 


1 




c(l + 


A) 


1 




c(l + 


A) 



lw">(0)]v(x T )- — [w'(0)]v'(x-), 



a 



[ w "'(0)}v(-x T ) + T ^-[w'(0)}v'(-x+) 

Using the connection function w(x) defined in l|1.3|l . we have 

K(0)] = w'(0+) - w'(0-) = 2(1 - aA) 

[w"(0)}=w"(0+)-w"(0-) = 

[w"'(0)] = w'"{0+) - w"'(0-) = 2(1 - a 3 A) 

These results lead directly to the following theorem. 

Theorem 4.1. The continuous eig enf unction v(x) on (—00,00) in has the 

following jumps in its first, second and third order derivatives at the boundary —x T 
and x T . 

(4.17) [v(x T )\ = 

(4.18) [v'(x T )j = 2a{ l~ a x A) v{x T ) 

(4.19) [v"{x T )\ = 2( "f~^ v(x T ) 

c(l + A) 

(4.20, = ^.(x,) + 

(4.21) W-x T )] = 

(4.22) [v'{-x T )\ = 2a ^- a x A \ {- XT ) 

(4.23) [v"(-x T )\ = zgMzil ^-^) 

(4.24) K(-x t)] = 2 -±^v(-x T ) - 2 *^v>(-xi). 
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ODE I ODE II ODE III 

v/x) v 2 (x) v 3 (x) 

I I I 

Xji 

Fig. 4.1. Valid ODEs on different sections and their solutions 

4.3. Eigenfunction symmetries. We define v%(x), V2(x) and ^(x) as the so- 
lutions of ODE I, ODE II and ODE III, respectively (see Fig. Hil l The three ODEs 
are all linear with constant coefficients. The continuous and bounded eigenfunction 
v(x) of l|2.11|l is defined as the following 

!vi(x), x € (— oo, — x T ], 

v 2 (x), x e [—x T , x T ], 
v s (x), x e [x T ,oo), 

and v\(x) matches v 2 {x) at — x T and V2(x) matches 113(2;) at x T . 
Lemma 4.2. The eigenfunction v(x) is either even or odd. 

Proof. By symmetry of ODE II, if V2{x) is a solution then v 2 (—x) is also a solution 

xj wu *, +■ v 2 (x) + v 2 {-x) v 2 (x)-v 2 {-x) 

Hence, both the even function and the odd function - 

2 2 

are solutions of ODE II. 
Let 

w(x - y)v 2 (y)dy. 

-XT 

If v 2 (x) is an even function, then since w(x) is even, T 2 {x) is also even. 

By the continuity of v(x) on 5ft, v(x T ) and v{—x T ) can be replaced by v 2 (x~) and 
v 2 (—x T ~), respectively. Thus the eigenvalue problem 12.11f) is 

(4.25) (1 + X)v(x) = w(x - Xt)^^± + w(x + x T )^^ + T 2 {x) 

c c 

Given that v 2 (x), w(x) and T 2 (x) are all even functions, from l|4.25|l . we see that v(x) 
is also even. Similar, we can show that v(x) is odd when v 2 (x) is odd. □ 

Lemma 4.3. The matching conditions at —x T are identical to those at x T when 
v{x) is an odd or an even function. 

Proof. This is shown with a direct calculation of the matching conditions of v'(x), 
v"(x) and v"'(x) at both —x T and x T . 

If v(x) is even, i.e. v(—x T ) = v(x T ) and v'(— x^) = —v'(x~), then defining the 
jump of v at x as [w(a;)] = v(x + ) — v(x~), the follow equalities are derived 

(4.26) [«(-it T )] = - [v(x T )} 

(4.27) [v'(-Xt)} = [v'(xt)} 

(4.28) W'{-Xt)\ - - [v"{xt)\ 

(4.29) W"{-Xt)\ = [v"'(xt)] 

Given the equalities (|4.26(l - (|4.29(l . a direct calculation shows that the matching con- 
ditions H4.21H4.24[1 at — Xt are equivalent to the matching conditions H4.17H4.20[1 
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When v(x) is odd, using the same approach, we can also justify that the matching 
conditions at —x T and x T are the same. □ 

4.4. ODE Solutions. ODEs I, II, and III are linear with constant coefficients 
and can be readily solved in terms of the parameters A, a, a, and u T . The eigenvalue 
A is specified when the solutions of the three ODEs are matched across the boundaries 
at x — —xt and x — xt- Solutions of ODE I are related to ODE III by a reflection 
x — > —x. By Lemma |4.3I the matching conditions at —x T are the same as those at 
x T . Thus matching solutions V2(x) of ODE II with solutions 113(2;) of ODE III across 
xt are sufficient to specify the eigenvalues of (|2.11() . The solution of ODE III is 

v 3 (x) = c 5 e~ ax + c 6 e~ x , 

where C5 and cq are constants. Notice that v 3 (x) exponentially decays to zero as 
x — ► 00, in accordance with the assumed properties of v(x). 

The solutions of ODE II will depend nontrivially on the parameters A, a, and a. 
The characteristic equation of ODE II is 

(1 + A)w 4 - Buj 2 + C = 0, 

where 

(4.30) B = {l + \)(a 2 + l) + 2a(l-aA) 
and 

(4.31) C= (1 + A)a 2 -2aa(A-a). 
The characteristic values are 



(4.32) uj 2 = 

where 



B± s/E 
2(1 + A) 



(4.33) A = B 2 - 4(1 + X)C 

= (a 2 - 1) 2 A 2 + 2(a 2 - 1) (a 2 - 1 - 2aAa - 2a)\ - 
(a 2 - 1)(1 - a 2 + 4a + AaAa) + 4a 2 (l - aA) 2 

Let Xb be the zero of B. If A is negative, l|4.32() shows that ODE II will have 
complex characteristic values. If A is positive, combinations of B and A yield either 
real or imaginary values. For fixed A, a and a, A is a parabola with a left zero A/ 
and a right zero A r . By Lemma Ij8.9|l and Lemma H8.10|l in the Appendix, cither 
A; < Xb < A,- and does not intersect with either branch of \fK or < A; and 
intersects with the left branch of VA. Table 14.11 and 14.21 describe all the possible 
structures of the characteristic values ±u>i and ±W2- There are three possible forms 
of solution V2(x): 1) both oj\ and u>2 are real; 2) both U)\ and L02 are complex; 3) u>\ 
is real and u>2 is imaginary. 

We notate the even solutions of ODE II as v^x) and the odd solutions as v^ix). 
When X > X r or Xj < X < Xi, both u}\ and u>2 are real. Thus 

(4.34) vl{x) = c 3 m(x) + Ci tM^M 

OJi — LO2 
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Table 4.1 

Characteristic value chart when A; < \g < A r 



1 


2 


3 


4 


5 


— 1 < A < A; 


A = A; 


A ; < A < X r 


A Ay- 


A > A r 


B < 


B < 


B > or B < 


s > 


B > 


A > 0, 


A = 


A < 


A = 


A > 


\B\ < VA 










lui real 


cui/2 imaginary 


oji,2 complex 


0^1,2 real 


0)1,2 real 


u>2 imaginary 




Wl = W2 


Wl = <^2 





Table 4.2 

Characteristic value chart when Xg < A; < A r 



1 


2 


3 


4 


5 


6 


-1< A < A 7 


A/ < A < A; 


A = A, 


A; < A < A r 


A Ar 


A > A r 


B < or 
S > 


B > 


B > 


B < 


B > 


B > 


A > 0, 
\B\ < VA 


A > 0, 

|s| > Va 


A = 


A < 


A = 


A > 


lji real 
L02 imaginary 


uj\,2 real 


lji : 2 real 

Wl = 0>2 


wi,2 complex 

Wl = U>2 


wi,2 real 

Oil = Ul2 


071,2 real 



where /zi(x) = e UlX + e^ lX , and ju 2 (ar) = e" 21 + e~ W2a: . We use because it is 

more convenient to resolve the degenerate case of u>i = lo 2 . As A — * A~, jJbx — ► /12 and 
e = — — * 0, l|4.34[) becomes 



Replacing e ca: by 1 + ex, e ex by 1 — ex and taking the limit as e — * yields 
(4.35) 



v e 2 (x) = c 3 (e Wl2; + e-^ x ) + c iX {e^ x - e"" 1 *) 
= 2C3 cosh pa; + 2c±x sinhpa; 



(|4.34|) approaches (|4.35|) as A — » A,7. It matches the solution v 2 (x) as A — * A+, which 
is given in l|4.37[) . 

Similarly, v 2 (x) can be written as 



Ul\ — UJ2 



When A/ < A < A r , u)\ and u>2 are complex. Let u)\ = p + iq, lu 2 = P — iq- When 
V2(x) is even, write v 2 (x) as 



(4.36) 



B . smqx 
v 2 (x) = ZC3 cos qx cosh px + ZC4 smhpa; 



As A — > A^ or A r , q — > 0, 

(4.37) ^C 2 -) * coshpx + 2C42; sinhpx 
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v 2 (x) can be written as 



2c3 cos qx smhpx 



sin qx 
2c4 coshpx 



where p — 
When 



2(1 + A) 

-1 < A < 



cos 



p : 



2(1 + A) 

u>i is real and w 2 is imaginary. 



■ sin 9 and — — arctan ■ 
2 

Let LL>2 - 



B 

iq, where 



2(1 + A)' 



then 



(4.38) 
(4.39) 



v 2 {x) = c 3 (e" i:E + e- UlS ) + 2c 4 cos(qx) 



v° 2 (x)=c 3 (e^* 



c ) + 2c. 



sin(ga;) 



Q 



5. Stability criteria. By theorem l4.1l vi (x) and vn(x) must match at — x^. and 
1*2(2;) and 1*3(2;) must match at x T . By property 14. 31 the matching conditions at — x T 
are same as the matching conditions at x T for v(x) even or odd. Therefore, it suffices 
to apply the matching condition to 1*2(0;) and 1*3(2;) at x T for the even and odd cases 
separately. This reduces the dimensionality of the resulting eigenvalue problem by a 
factor of two. In general, the matching conditions can be written as 



Tl 



\v{x T )} = 1*3(2;+) - v 2 (x T ) 
[v'(x T )] = v' 3 (x+) - v' 2 (x-) 

[v"(x T )} =v%(x+)-v%(x-) 
[v'"(x T )\ = v'i'(x+) ~ v' 2 "(x- 





2a(l 



aA) 



1 + A 
2{aA-l] 



v{x T ) 



c(l + A) 
2(1 - a 3 A) 

c(l + A) 



v(x T ) 



v{xi 



2a(aA- 1) 

TTa 



v'(x-) 



where v(x T ) — ^3(0;+) and v'(x~) = v' 2 (x~). 

A given stationary pulse solution uq(x) will be specified by a set of parameters 
a, A, a, x T , and u T . The eigenvalues A that determine stability of pulse solutions 
are given by system Tl. To compute these eigenvalues, we require the appropriate 
form of the eigcnfunctions v 2 (x) and ^3(2;). We do so by finding characteristic values 
(|4.32|) corresponding to the parameters specifying the given stationary pulse solution. 
We expediate this process by calculating the constants B l|4.30|l and C (|4.31|) . then 
using Table 14.11 and Table 14.21 to deduce the characteristic value types. We then 
substitute the appropriate form for ^2(2;) and 1*3(2;) into T\ where coefficients C3 and 
C4 in t>2(2; T ) and C5 and cq in V3(x T ) are unknown. We replace v(x T ) by 1*3(2;^) and 
v'(x~) by v' 2 (ly). This results in a 4 x 4 homogeneous linear system with 4 unknown 
free parameters C3, C4, C5, cq. We must do this for both even and odd eigenfunctions 
resulting in two separate linear systems that must be solved. 

The coefficient matrix of this system must be singular for a non-trivial solution 
(03, C4, C5, cq). Hence, the determinant D(X) of the coefficient matrix must be zero. 
Thus, the solution of D(X) =0 is an eigenvalue and it determines the stability of the 
stationary solution. If there exists a A such that < A < A& and D(X) = 0, then 
the standing pulse is unstable. If there is no positive A such that < A < A& and 
D(A) — 0, the standing pulse is stable. Our determinant D(X) for stability is similar 
to the Evans Function [131 ITSl IT71 ITS] . 
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5.1. Stability of the small and large pulse. Two single-pulse solutions were 
shown to exist in the accompanying paper |32J for parameters a = 2.4, A = 2.8, 
a = 0.22, u T — 0.400273 and = 1. The large pulse has a higher amplitude and 
larger width and is denoted by u (x). The small pulse is slightly above threshold and 
much narrower than v}(x) and is denoted by u s (x). The explicit forms are given by 

, J 0.665 cos(0.31a;) cosh(1.49x) - 3.78 sin(0.31a;) smh(1.49x) + 0.33, x g [~x T ,x T ] 
U ^ X '~ \ 6.237e" 2 ' 4W - 1.604e _w , otherwise 

where x T = 0.683035. 

0.22 cos(0.31x) cosh(1.49s) - 8.03 sin(0.31a;) sinh(1.49a;) + 0.33, x € [-x T ,x T ] 
1.203e _2 4|a:| - 0.416e _w , otherwise 

where x T = 0.202447. 

We first calculate the upper bound for the eigenvalue Af, which is different for the 
large pulse and small pulse because Af, depends on x T . Let X\ be the upper bound for 
the large pulse and A^ be the upper bound for the small pulse. For the parameter set 
a = 2.4, A = 2.8, a = 0.22, u T = 0.400273, the upper bounds are Aj, = 1.25917 and 
A* = 1.66628. 

For the above set of parameters, v^(x) always has the following form. 

v 3 (x) = c 5 e~ ax + c 6 e~ x . 

The form of Ua(a;) depends on uj\ and u>i. For this specific set of parameters, the 
left and right solutions of A lj4~3"3"|) are A; = -0.627692 and A r = 0.192861. When 
< A < X r , both u>\ and u>2 are complex, implying 



V2{x) 



. . smqx ... 

v 2 (x) = 2c3 cosqx cosh pa; -I- 2c4 smhpx V2\x) is even 

Voix) — 2c3 cos qx sinh px — 2c4 cosh pec v-iix) is odd 

q 



where p, q are real, and C3, C4 are unknown. 

Substituting (v 2 (x)) and Vq(x) into system Tl, results in an unwieldy 4x4 

linear system in C3, C4, C5 and cq. We use Mathematica |55] to calculate the determi- 
nant of the coefficient matrix as a function of A. 

When 0.192861 = A r < A < AJ, = 1.25917, wi, 2 is real, and v 2 (x) has the form 



v 2 (x) 



c^ lX + e UJlX ) + c 4 - v 2 {x) is even 

LUl - LJ2 



c 3 (e" i:E - e-^ x ) - c 4 - i : w 2 (x)isodd 

W\ - U2 



Figure I5~T1 gives a plot of D(X) on the domain [0, Af,], combining the regimes where 
wi,2 is real and complex. We see that there is no positive A that satisfies D(\) = 0. 
Figure l5~2l shows D(X) for odd v(x). We see that D(\) = only when A = 0, which 
is consistent with Theorem 18. 61 The lack of a positive zero of D(X) indicates that the 
large pulse is stable. 

For the same set of parameters, {a = 2.4, A = 2.8, a = 0.22, u T = 0.400273}, the 
upper bound of the small pulse is Ajj = 1.66628 Repeating the same procedure as 
for the large pulse, we plot D{\) for both v^ix) and (Figs 15.31 and l5.4|l . The 

existence of a positive eigenvalue A = A* satisfying D(X*) = in Fig. 15.31 implies the 
instability of the small single-pulse. The plot of D(A) corresponding to v°(x) in figure 
15.41 identifies the zero eigenvalue. 
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D(X) 




Fig. 5.1. Plot of D(X) for large single-pulse v}(x) when Vzix) 



is even, a 



2.4, A = 2.; 



a = 0.22, u T = 0.400273, x T = 
A such that D(X) =0, A < \\- 



0.683035, A r = 0.192861, \[ = 1.25917. There is no positive 




5.2. Stability and instability of for different gain a. For both the large 
single-pulses and small single-pulses, D(X) is monotonically increasing (See Fig 15.51 
and 15.61) . However, D(0) for small pulses is negative. As A increases, D(\) crosses 
the A-axis and becomes positive. Therefore, -D(A) has a positive zero. For the large 
pulse, -D(O) is positive and D(X) has no positive zero. We follow D(0) for a range of 
a E (0.22, 0.59) in Fig 15. 71 and find that D(0) is always negative for small pulses and 
positive for large pulses. Hence, the large pulses are stable and the small pulses are 
unstable in this range. 
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Fig. 5.4. Plot of D(X) for small single-pulse vf(x) when V2{x) is odd. a — 2.4, A = 2.8, 
a = 0.22, u T = 0.400273, x T = 0.683035, A r = 0.192861, X B b = 1.66628. There is no positive 
X such that D(X) — 0, A < A^. When V2(x) is odd, D(X) = at X — identifies the zero 
eigenvalue. 



5.3. Stability of the dimple-pulse u d (x) and the instability of the third 
pulse. When there are only two single-pulses, the large pulse could be a dimple-pulse 
instead of a single-pulse. This dimple pulse has the same stability properties as a large 
pulse. The parameter set a = 2.4, A = 2.8, a = 0.22, u T = 0.18, and x T = 2.048246 
corresponds to a dimple pulse. Carrying out the stability calculation yields D(\) 
shown in Figs. 15.81 and 15.91 We see that there is no zero crossing and thus the dimple 
pulse is stable. This is true for all dimple pulses we tested in this category. 

As shown in |24) and |25 |. for certain parameter regimes, there can be more than 
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D(X) 



50 




0.5 1 



Fig. 5.5. Plots of D(X) for large single-pulses with different gain a. a — 2.4, A = 2.. 
a = 0.22, ur = 0.400273. 



D(X) 



40 



-40 



-80 




Fig. 5.6. Plots of D(X) for small single- 
a = 0.22, u T = 0.400273. 



es with different gain a. a = 2.4, A — 2.8, 



two coexisting pulses. When there are three pulses, the third pulse can be either 
a single-pulse or a dimple-pulse. For example, when A = 2.8, a = 2.2, a = 0.8, 
u T = 0.2, the third pulse is the single-pulse 

Jl.28 cos(0.47a;) cosh(1.2a;) + 1.27 sin(0.47:r) sinh(1.2a;) + 0.8129, x G \-x T , x T ] 
U ^ X ' [I98.78e 2|a;| - 15.15e- |a;| , otherwise 

where x T — 2.20629. D{\) shown in Fig. 15.101 indicates that this pulse is unstable. 
When a = 2.6, A = 2.8, a — 0.6178, u T — 0.063, the third pulse is the dimple pulse 



u(x)= 



0.35 cos(l. 112a;) cosh(1.112x) + 0.24 sin(1.112x) sinh(l. 112a;) + 0.163, x G \-x T ,x T ] 
232.89e 26|a:| -9.31e~ w , otherwise 
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0.22 



0.59 



-40 - 




Fig. 5.7. Plots of D(0) for both large single-pulses (blue branch) and small single-pulse 
(red branch) with a G (0.22,0.59). a = 2.4, A = 2.8., u T = 0.400273. 

D(k) 




\ 1 \ 



Fig. 5.8. Plot of D(X) for dimple-pulse when Vi{x) is even a — 2.4, A — 2.8, a — 0.22, 
x T = 2.048246, X r = 0.192861, = 2.48147. There is no positive X such that D(X) = 0. 



where x T = 1.98232. As seen in Fig. 15. Ill D(X) crosses zero for a positive A indicating 
that it is unstable. In all the cases that we have examined, we find that the third 
pulse is unstable. 

6. Double-pulse and its stability. For certain parameter regimes, there can 
be double-pulse solutions which have two disjoint open and finite intervals for which 
the synaptic input u(x) is above threshold |241 1251 132| . An example is shown in 
Fig. 16. 1|) . We consider symmetric double-pulses that satisfy the equation 



/Xx r-XQ, 
w(x-y)f[u(y)]dy+ J w(x - y)f[u{y)]dy 
-X2 J X\ 
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Fig. 5.9. Plot of D(X) for dimple-pulse when V2(x) is odd. a — 2.4, A — 2.8, a — 0.22, 
u T = 0.18, x T = 2.048246, A r = 0.192861, A^ = 2.48147. There is no positive A such 
that D(X) = 0, A < A^ . When Vzix) is odd, D(\) does identify the zero eigenvalue because 
-D(A) =0 at A = 0. This is consistent with theorem. 




Fig. 5.10. Plot of D(\) for the third pulse (a single-pulse) when V2(x) is even, a — 2.2, 
A = 2.8, a = 0.8, u T = 0.2, x T = 2.0629, c = 2.75017, D(0) = -0.0153. There is a positive 
A such that D(X) — 0. 



where x\_i > 0. Thus u > u T for x G (xx,X2) U (— x%,— x\), u = u T for x — 
—X2, —x\,xi,X2, and u < u T outside of these regions and approaches zero as x — > oo. 
We show their existence in |23] and |25j . 

Linearizing the dynamical neural field equation (|1.1|) around a stationary double- 
pulse solution u{x) gives eigenvalue equation: 

/„ «s / ,x / x / -,v(xi) , . v(— Xi) , ^v(X2) 
(6.2) (1 + X)v(x) = w(x - xi)-^- + w(x + xi)— + w(x - x 2 )^- JlL 

Ci Cl c 2 




Fig. 5.11. Plot of D(\) for the third pulse (a dimple-pulse) when V2(x) is even, a — 2.6, 
A = 2.8, a = 0.6178, u T = 0.063, x T = 1.98232, c = 2.21523, D(0) = -0.094. There is a 
positive A such that D(X) — 0. 

LI 

0.8 r 




-0.3 



Fig. 6.1. Double-pulse for Amari case in which a = 0. A = 2.8, a = 2.6, a = 0, 
Ut = 0.26, aci = 0.279525, x 2 = 1.20521. 



+ w(j: + a;2)— — + a ( / w(x — y)v(y)dy + I w(x — y)v{y)dy 



C 2 

The eigenvalues A of (|6.2|l possess the same properties as those of the eigenvalue 
equation for the single-pulse solutions. 

For simplicity, we consider the Amari case in which a = 0. The solution of (|6.2|l 
for a > would involve a long calculation. For a = 0, the eigenvalue equation (|6.2|l 
becomes 

(6.3) (1 + X)v(x) = w(x — x\) - 1 ^ + mix + xi)— — 

Ci Ci 

w(x — X2) h iv[x + x 2 )- 



C2 C 2 
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where C\ = v!(x\) and c 2 = u'(— x^)- Then v!{— X\) — —ci and u'ixz) = — &2- Using 
an approach similar to Theorem 18. 41 in the Appendix, we can show that A is real. By 
taking the derivative of jSJ , we can also show that zero is an eigenvalue of system 
(|6.3(l . and the corresponding eigenfunction is u'(x). 

Setting x = X\, x = — x\, x — x-i and x = —x-i in (|6.3|) gives a 4-dimensional 
system 



/ w(0) - A - 1 

Cl 

w(2xi) 



(6.4) 



w(2xi) w(xi — X2) w(xi + X2) \ 



Cl 



10 (0) 



A 10(2:1 + £2) w(a;i — x%) 



Cl 



Cl 



f'2 



C2 



w(a;i — K2) to(a?i + x?) w(0) — 1 — A 



w(2a;2) 



Cl 



Cl 



C2 



C2 



V 



w(xi + a: 2 ) 

Cl 



w(xi 



X2 



Cl 



w(2gg) 

C2 



io (0) 



/ ^1) \ 

v(—xi) 
v(x 2 ) 



= 



C2 



The determinant -D(A) of coefficient matrix in system (16.411 is a fourth order poly- 
nomial. Since zero is an eigenvalue, then D(X) = Xd(X), where d(X) is a third order 
polynomial. Consequently, the stability of the stationary solution u(x) is determined 
by the roots of a third order polynomial d(X), which can be found numerically. We 
computed d(X) for the two double pulses shown in Fig. 16.11 Figure IFT21 shows a plot 
of the third order polynomial d(X) for the small double-pulse. It has three positive 
zeros indicating instability. The plot of d(X) for the large double-pulse as shown in 
Fig. l6.3fl has two positive zeros. Therefore, both the small and large double-pulses are 
unstable. We have not found any stable double-pulses for any parameter sets that we 
tested. However, we have not fully investigated the parameter space of A, a and u T . 

d(k) 




Fig. 6.2. Plot of polynomial d(A) for the small double-pulse shown in Fig. 16'. 1\ 



7. Discussion. Our results show that although many types of pulse-solutions 
are possible, only the family of large pulses and associated dimple pulses are stable. 
For the situation of three coexisting pulses, the third and largest pulse is always 
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Fig. 6.3. Plot of polynomial d(X) for the large double-pulse shown in Fia \6.1\ 



unstable. It is possible that more than three pulses can coexist although we did not 
investigate situations beyond three. The double pulses we examined were not stable 
in accordance with previous work [32] • 

The caveat is that we were only able to examine specific examples individually 
or over limited parameter ranges. Although we have an analytical expression for the 
eigenvalues the length of these expressions makes them difficult to analyze. As a 
result, we were unable to make as strong a claim as Amari who showed that large 
pulses are always stable and small pulses are always unstable j3J- It may be possible 
to find some patterns in the expressions to make more general deductions. From our 
parameter explorations, we were unable to find stable pulse solutions other than the 
large and associated dimple pulse. 

We wish to note that numerical simulations on discretized lattices can give mis- 
leading results regarding the stability and existence of pulse solutions of the associated 
continuum neural field equation. We conducted some numerical experiments using a 
discretization of the neural field equation l|l.lfl and to our surprise we were able to 
easily find examples of stable dimple and double pulses even though the continuum 
analogue shows that these solutions either do not exist or cannot be stable. The 
resolution to this paradox is that a discrete lattice may stabilize solutions that are 
marginally stable in the continuum case. 

Consider the Amari neural network equation consisting of N neurons 



where w(i — j) is given by (|1.3fl . O(-) is the Heaviside function, and Aa; gives the 
discretization mesh size. For an initial condition for which Uj > u T on a contiguous 
set of points {i...k} and k — i is less than the expected width of the large pulse in the 
analogous continuum neural field equation, the numerical solution converges towards 
the expected large pulse solution. However, if the initial set of points is larger than 
the width of the large pulse (we have not fully investigated how much larger it needs 
to be), then there is a possibility that the simulation will converge towards an entirely 
different state. 




N 



(7.1) 



3=0 
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u(x) 




0.2 - 

u T - 




Fig. 7.1. Result of numerical simulation of \7. If for parameters N = 200, Ai = 0.1, A = 1.8, 
a = 1.6, and u T = 0.124. The arbitrary discretization length scale is chosen so that x = O.li. 



For example, a numerical simulation of the parameter set N — 200, Ax = 0.1, 
A = 1.8, a = 1.6, and u T — 0.124 with an initial condition ui = 1 for i E 50. ..150, 
converges to a stable dimple-pulse state shown in Fig. 17.11 Different initial domains 
will lead to different attracting states where the width is close to the initial domain 
width. For a large enough initial domain, the dimple pulse will break into a stable 
double-pulse. Increasing the initial domain can lead to increasingly higher number 
stable multiple pulses. 

We can show that these states do not exist in the analogous continuum neural 
field equation Consider a stationary pulse solution of (|1.1|) for a = 0. A pulse of width 
x T satisfies 

(7.2) u(x) = 4>(x,x T ), 
where 

(7.3) ct>(x,x T ) = f T Ae~ a \ x - y \ - e~\ x - y \dy. 

J —XT 

The pulse can exist if it satisfies the existence condition 

(7.4) u T = 4>{x T , x T ) 

from which the width be obtained. A plot of the existence condition is shown 

in Fig. O 

It is immediately apparent that the large pulse does not exist. The existence 
function approaches u = u T from above for large enough x T . While it is very close 
to u T it never crosses it. However, for the analogous discretized equation H7.1|) . the 
discrete mesh can break the symmetry of this nearly marginal mode and result in a 
family of stable pulse solutions for arbitrary widths larger than a given width. 

This effect can be intuitively understood by examining Fig. 17.11 The neurons 
immediately adjacent to the edge of the pulse are significantly below threshold and 
thus have no effect on the rest of the network. A perturbation on the order of the 
distance they are below threshold would be necessary to cause these neurons to fire 
and influence the network. In the continuum equation, the neurons on the boundary of 
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Fig. 7.2. Existence condition for pulse solutions of neural field equation |7.£| ) for parameters 
A = 1.8, a = 1.6, and u T = 0.124. 



the pulse are precisely at threshold. Arbitrarily small perturbations can push the field 
above threshold and influence the other neurons. A stable pulse must withstand these 
edge perturbations. Discretization eliminates these destabilizing edge perturbation 
effects. 

We can make a simple estimate of how fine the discretization mesh must be in 
order for these discrete affects to disappear. The distance the neuron adjacent to the 

pulse is below threshold is approximately given by d x (j)( X X rp , X T )dx ~ (A - l)dx. 

For the parameter set of our simulation, the continuum existence condition shows 
that 4>{x T ,—x T ) — u T > 0.001. Thus to eliminate the discreteness effect, we require 
the adjacent neuron to be above threshold, i.e. (A — l)dx < 0.001 as it would be in 
the continuum case. This leads to an estimate of dx < 0.00125. Hence, for a domain 
of dimension x > 20, a network size of N > 16, 000 is necessary to eliminate the 
discreteness effect. 

Biological neural networks are inherently discrete. Thus this discreteness effect 
may be exploited by the brain to stabilize localized excitations. Our numerical simu- 
lation is an example of a discretized line attractor |51| where the width of the pulse 
is determined by the initial condition. Although, the discrete network may have a 
richer structure, this does not imply that the study of continuum neural field equa- 
tions are not necessary. Field equations lend themselves more readily to analysis and 
many insights into the structure and properties of neural networks have been gained 
by studying them. We suggest that studies combining neural field equations, discrete 
neural network equations and biophysically based spiking neurons may be a fruitful 
way to uncover the dynamics of these systems. 

8. Appendix 1: Properties of the eigenvalue problem. We prove some 
properties of the the eigenvalue problem l|2.11|l with the connection function given by 
(|1.3(l . First consider functions 

<h (*) = 7^ [°° e ~ alx ~ vl ( F u + ®u)v{y)dy 
2a J-oo 
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fa{x) = - e~^y\(F u + e u )v(y)dy 

^ J -oo 

where F(u) = a(u — u T ), 0(u) is the Heaviside function, and subscript denotes partial 
differentiation. 

Lemma 8.f . The eigenfunction v(x) satisfies 

(1 + X)v = 2(aAfa - fa) 



Proof. 

„ ., , ,v(x T ) .v(—x T ) f XT . . . . , 

(1 + X)v — w(x — x T ) \-w(x + x T ) \- a / w(x — y)v(y)dy 

c c J_ XT 

5(x — x T ) + S(x + x T ) f°° 
w{x-y) v(y)dy+ w(x - y)F u v(y)dy 



(x - y)Q u v(y)dy + / w(x - y)F u v(y)dy 



= A e- a ^(F u +e u )v(y)dy- \ e~^{F u + B u )v(y)dy 

J — OO 

= 2{aAfa - fa) 



□ 



Lemma 8.2. Functions fa and fa satisfy 
(8.1) -4>'{ + a 2 fa = {F u + e u )v 



(8.2) 



-fai + a 2 fa = (F u + e u )v 



Proof. The second order derivative of (x) is 

r px poo 

(8.3) fa = - / e- a ^- y \F u +Q u )vdy+ e a ^(F u + <3 u )vdy 

* LJ — oo Jx 

- (f u + e u )v 

-(H3J + a? fa (x) yields 

(8.4) - 4l( + a 2 fa = (F u + Q u )v. 

—<j>2 + a 2 fa — (F u + Q u )v can be obtained in the same fashion. □ 

Lemma 8.3. lim fa % = and lim fa 2 = provided that v(x) is bounded on 

x — >±oo ' x — >±oo ' 

(— oo. oo) and exponentially decays to zero as x — > ±oo 
Proof. When x » x T 



fa(x) 



1 

2a 



ne ax I e ay v(y)dy + e 

~X T 



-a(x-x T ) v{Xt) -a(x+x T ) v (~ x t) 



Hence, lim fa = provided that v(x) is bounded on [— x T ,x T ]. 

x — >oo 

When x << —x T < 0, as x — » — oo 

fXT 



fa(x) 



1 

2a 



e- ay v(y)dy + e ^-T) v _^A + ga^^jHLEl) 

c c 
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As x — > oo, 



e -(»-») (F u + O u )vdy + / e «(»-^ (F u + Q u )vdy 



lim 0; = lim (-i / 



lim 

a: — *-oo 



a 



e ay dy + e 



As x — > —00 



lim cf>[ = _ lim _ ( - i I e a{x - y \F u + Q u )vdy 



x — > — 00 x — > — 00 



= lim < 

2 



= 



Similarly, one can prove that lim </>2 = and lim (f>' 2 = 0. Therefore, 

X — >±OC X — >±OG 

lim = and lim (/)[ 2 = 0. □ 

x — >±oo ' x — >±oo ' 

Theorem 8.4. The eigenvalue A in \2.11\l is always real. 
Proof. Using the results of Lemma T8.2I aAcbi (18.11) — 60 l|8.2|l gives 

(8.5) aAM-<t>" + a 2 4>i) ~ M~4>2 + <h) = {F u + G u )v(aAfa - fa) 
where <pi^2 are the complex conjugates of <p±.2- Integration by parts gives 



b'(dx 



\§\dx 



\4{fdx 

-00 



and similarly 



b' 2 f dx. ^,From Lemma l8~T1 



-(l + X)v = aA4> 1 - ( j) 2 
~(l + X)v = aA(pi-<p2 



Integrating both sides of l|8.5|l gives 



(8.6) 



Using 



aA 



i>'i\ dx + a 2 I \<j>i\ 2 dx 



\v\ <d u dx 



• 2 \'dx+ I \<f> 2 \' dx ) = -(1 + \) I \v\ 2 (F u + e u )dx 



1 f°° 1 / 

- / lul (S(x ~ x T ) + 5(x + x T j) dx — — \\v(x t )\ + \v(— x T )V 

c . „ c V 
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in H8.7fl and rearranging gives 



(8.7) 2( 1 + A ) = 



F u \v\ 2 dx + i (K* T )| 2 + |„(-^)1 2 ) 



The right-hand side of (|8.7(l is real, therefore A is real. □ 

Theorem 8.5. The eigenvalue X in ^2.11}) is bounded by X b = _|_ 2ak 1 x T — 1 
where fco is the maximum of \w(x)\ on [0, 2x T ] and \w(x — y)\ < k\ for all (x,y) G 
J x J, where J — [— x T , x T ] . 

Proof. We write the eigenvalue problem H2.11J) as 

(8.8) (l + X)v = Lv 

where operator L is defined as (|2.12|) . 

Function w(x — y) is continuous on square J x J. We take the norm of both sides 

of 



(1 + A)H = ||HI 



with norm 



max 

x6J 



Thus 



IMI 



umx — x T ) \-w(x + x T ) \- a / w(x — y)v(y)ay 

c c 



< max 

x£j 



w(x — X T ) 



v(x T ) 



max 

xGJ 



w(x + X T ) 



J -X T 

v{—x T ) 



max 

xeJ 



a / w(x - y)v(y)dy 



< \w(x — x T )\ h \w{x + x T )\ ha|ju|| / max \w(x — y)\ dy 



x£j 



\\v(x)\\ 

< 2fc + 2ak 1 x T \\v(x) 



where 



fco = max \w(x — x T )\ = max \w(x + x T )\ 

xGJ x£j 

since w(x) is symmetric and \w(x — y)\ < k\ for all (x, y) G J x J. Therefore 



(1 + A)||u(a;)|| = ||Li;(a;)|| < 2fc 



|u(x)|| 



2afcix T |ju(a;)|| 



leading to 



A < — - + 2afcia; T — 1 = A/,, 
c 



Theorem 8.6. X — is an eigenvalue. 
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Proof. Consider the equilibrium equation 

/>oo 

u(x) = / w(x - y)f[u(y)]dy 



w(x-y){a[u(y)-u T ] + 1} dy 

-XT 

where u(x) is a stationary standing pulse solution. After a change of variables p 
x — y, l|8.9l) becomes 

i-X + Xt 

(8.10) u(x) = / w(p) {a [u(x - p) - u T ] + 1} dp 

J x — xt 

Differentiating l|8.10|) with respect to x yields 

(8.11) u'(x) = w(x + x T ) [a(u(— x T ) — u T ) + 1] — w(x — x T ) [a(u(x T ) — u T ) + 1] 

px+xt 

w(p)u' (x — p)dp 

1 x — xt 

Since u{— x T ) = u(x T )u T and u'(—x T ) = c = —u'(x T ), 

u(x)=w(x + x T ) w(x — x T ) ha / w{p)u (x — p)dp 

C C J X—XT 

u'(x ) u'(-x ) f XT 

(8.12) =w(x-x T ) T +w(x + x T )— —+a w(x - y)u' (y)dy 



^ ^ J -X T 

(|8.12() is the eigenvalue problem l|2.11|l with eigenvalue A satisfying 1 + A = 1, resulting 
in A = 0. The corresponding eigenfunction is u'{x). Therefore, A = is an eigenvalue 
of l|2.11l) corresponding to eigenfunction u'(x). □ 
Theorem 8.7. Consider the operator 

(8.13) L = T 1 +T 2 , 

where 

Ti(v(x)) — w(x — x T ) - T ^ + w(x + x T )— — , T\ : C \—x T , x T ] — > C [— x T , x T ] 

c c 

/XT 
w(x - y)v(y)dy, T 2 : C [—x T , x T ] — * C [~x T , x T ] 

-XT 

Both T\ and T 2 and hence L are compact operators. 

Proof. It is obvious that both T\ and T 2 are linear operators. The boundedness 
of Ti follow from 



\T\v\\ = max 



v(x T ) M-x T ) 
w(x — x T ) V w(x + X T ) 



c c 



< \w(x — Xt)\ h \w[x + Xt)\ 

c c 

< 2fco M 

c 

Let v n be any bounded sequence in C [—x T , x T ] and \\v n \\ < cq for all n. Let y\ — T\v n . 
Then ||y^|j < ||Ti||||u„||. Hence sequence y\ is bounded. Since w(x,t) = w(x — t) is 
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continuous on J x J and J x J is compact, w is uniformly continuous on J x J. Hence, 
for any given t\ > 0, there is a Si > such that for t = x T and all X\,x% € J satisfying 
\x\ - x 2 \<8 1 

Q 

\w(xi - x T ) - w(x 2 - x T )\ < — ei. 

2c 

Consequently, for x%, x 2 as before and every n, one can obtain 

v„(x T ) 



IVnixi) ~ Vn{x 2 )\ = 



[w(xi — X T ) — Ul(x 2 — X T ) 



Vji\ x T j 



+ [w(xi + X T ) — W(X2 + X 7 



< \w(xi - x T ) - \w(x 2 - x T )\ — + \w(xi + x T ) - w(x 2 + x T )\ — 

c c 

c cq c c 
2cq c 2cq c 



Boundedness of T 7 follows from 



1^2^11 < H^ll max 



\w(x - t)\dt 



Similarly, let y\ — T 2 v n . Then is bounded. For any given e 2 > 0, there is a S 2 > 
such that for any t € J and all x\,x 2 G J satisfying \x\ — x 2 \ < S 2 

c 2 

\w(xi -t)- w(x 2 - t)\ < 



\yl(xi) - yl(x 2 )\ 



[w(x\ — t) — w(x 2 — t)]v n (t)dt 



< 2x 7 



(2 



2x t cq 



(2 



This proves the equicontinuity of {y„} and {y„}- By Ascoli's theorem, both se- 
quences have convergent subsequences. v n is an arbitrary bounded sequence and 
y\ = T\v ril y^ — T 2 v n . The compactness of T\ and T 2 follows from the criterion that 
an operator is compact if and only if it maps every bounded sequence x n in X onto 
a sequence Tx n in Y which has a convergent subsequence. □ 

Theorem 8.8. A = — 1 is the only possible accumulation point of the eigenvalues 
of L and every spectral value A —1 of L is an eigenvalue of L. Thus the only 
possible essential spectrum of compact operator L is at A = — 1. 

Proof. Let 7 — (1 + A), the eigenvalue problem becomes 

jv(x) = Lv(x), 

and the linear operator L is compact on the normed space C[— x T ,x T ] . 7 is the 
eigenvalue of operator L. The spectrum of a compact operator is a countable set with 
no accumulation point different from zero. Each nonzero member of the spectrum 
is an eigenvalue of the compact operator with finite multiplicity [JHIEHI- Therefore, 
the only possible point of accumulation for the spectrum set of compact operator L 
is 7 = 0, i.e., A = —1 and every spectral value A 7^ —1 of L is an eigenvalue of L. 
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This suggests that the only possible essential spectrum is at A = —1. All the spectral 
values A such that A > — 1 are eigenvalues. □ 

Lemma 8.9. The zero of B, X B , obeys -1 < X B < A r . For the case a 3 > A, 
Xi < Xb < X r , and for the case a 3 < A, Xb < A; < A r . 

Proof. Set 



B = (1 + A)(a 2 + 1) + 2a(l - aA) = 



The zero of B is 



a 2 + 1 + 2a - 2aAa , 2a(aA - 1) 

Xb = m = ~ l + 2,i > ~ L 

a 2 + 1 er 2 + 1 

A is a quadratic function in A and it has two zeros. The left zero is 

1 - a 2 + 2aAa + 2a - 4a\/aA 



A, = 



a 2 - 1 



The right zero is 



_ 1 - a 2 + 2aAa + 2a + AayfaA 
A r — 



a 2 - 1 

The difference between A r and A# is 



. . 4aa(a + A) + AaVaA(a 2 + 1) 

A r - X B = 2 — i > 

a 4 — 1 

Therefore —1 < Xb < X r . 

The difference between X B and A ; is X B — A ; = ^ a (^~®^ — j-)( a VaA) ^ 

a 4 — 1 

sign of Xb — A; depends on a 2 — \fa~A. If a 2 — \/aA is positive, i.e. o? > A, then 
Xi < Xb < X r . If a 2 — \J~aA is negative, i.e., a? < A, then As < A; < A r . □ 

Lemma 8.10. (i) For a 3 > A and Xi < Xb < X r , B does not intersect the left 
branch or the right branch of \/A. (ii) For a 3 < A and Xb < A; < A r , B intersects 
only the left branch of VK once at Xi . 

Proof. It is not difficult to see that B does not intersect the right branch of V~K 

for both (i) and (ii). \/A is linear in A with slope a 2 — 1 for large A. The slope of B 

is a 2 + 1. Both a 2 — 1 and a 2 + 1 are positive and a 2 + 1 > a 2 — 1, thus B and the 

right branch of ^fK never meet. When A; < X B < A r , B < for A < X B and \fK > 

for A < A; < Xb- Therefore B and y/A never intersect. In (ii), B intersects the left 

i i r nr x 2^4a — 2act — a 

branch of VA at A/ = . □ 

a 
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